A1 = readmatrix('10Vcool_30K.txt')
A2 = readmatrix('10Vcool_40K.txt')
A3 = readmatrix('10Vcool_50K.txt')
A4 = readmatrix('10Vcool_60K.txt')

x1 = A1(:,1)
y1 = A1(:,2)*1e12

x2 = A2(:,1)
y2 = A2(:,2)*1e12

x3 = A3(:,1)
y3 = A3(:,2)*1e12

x4 = A4(:,1)
y4 = A4(:,2)*1e12

hold on
xlim([0, 40])
plot(x1,y1,'Linewidth',1,'Color',[0.2 0 0.8])
plot(x2,y2,'Linewidth',1,'Color',[0.4 0 0.6])
plot(x3,y3,'Linewidth',1,'Color',[0.6 0 0.4])
plot(x4,y4,'Linewidth',1,'Color',[0.8 0 0.2])
xlabel('Polarization')
ylabel('Photocurrent [pA]')
ylim([-0.2,1.5])
yticks(0:0.4:1.2)
xticklabels({})
yticklabels({})
%set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
hold off

set(gcf,'Position',[0, 0, 200, 185])